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Abstract 

Based on the acoustic analog, we investigate both of the effects : disorder and (Pauli- 
blocking) interaction to the possible localization in electron gases by using the quantum 
discrete kinetic model. Effects of the disorder (or free-orientation : 9 which is related to the 
relative direction of scattering of particles w.r.t. to the normal of the propagating plane- 
wave front) which is introduced into the Uhling-Uhlenbeck equations together with those 
of the Pauli-blocking arc presented. We obtain the possible phase diagram (related to the 
strength of disorders and the mean free path) which resembles qualitatively that proposed 
by Abrahams 1 . 

PACS : 71.10.Hf; 71.30.+h; 73.43.Nq; 73.20.Jc; 73.20.Fz 

An intriguing aspect of quantum mechanics is that even at absolute zero temperature quan- 
tum fluctuations prevail in a system, whereas all thermal fluctuations are frozen out. These 
quantum fluctuations are able to induce a macroscopic phase transition in the ground state of 
a many-body system, when the ratio of two competing terms in the underlying Hamiltonian is 
varied across a critical value [1-4]. The instability of the quantum critical behavior with respect 
to (w.r.t.) the disorder can be interpreted as a signal for phenomena of localization. Recently 
some of the new developments in the localization problem [5] have become a major theme in 
the condense matter research. One example is the strongly-interacting electron (non-Fermi) liq- 
uid with different strengths of disorder [6]. Interesting issues are the quantum phase transition 
(QPT) and quantum critical point (QCP). Both effects of disorder and interaction are closely 
relevant to the weak and strong localization [8, 9]. They are then related to the metal-insulator 
transition in two-dimensions (2D). Although most of theories proposed before are based on the 
Fermi liquid behavior, new insights could be obtained considering the possible analogy with 
superconducting transition which might be related to the bosonic system [2-3]. 
In the last two decades a variety of metals have been discovered which display thermodynamic 
and transport properties at low temperatures which are fundamentally different from those of 
the usual metallic systems which are well described by the Landau Fermi-liquid theory. The re- 
sistivity in a variety of high mobility 2D electron/hole systems is seen experimentally to exhibit 
a number of interesting anomalies that don't as yet have an adequate theoretical understand- 
ing. Note that, at sufficiently low electron densities, an ideal two-dimensional electron systems 
becomes strongly correlated, because the kinetic energy is overpowered by energy of electron- 
electron interactions (exchange and correlation energy). The interaction strength is normally 
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described by the Wigner-Seitz radius, r s = 1/ (Ttn s ) l l 2 aB (where n s is the electron density and 
as is the effective Bohr radius in semiconductor). Till now, for rather large r s (S> 1) and to- 
gether with disorder, the nature of this metal-insulator transition (MIT) remains the subject of 
ongoing debate [6,8-9]. 

Motivated by the analogy between electrons in periodic or disordered metals and waves in clas- 
sical acoustical systems [10-12] an investigation for observing possible QCP in MIT or relevant 
localization [8] using the quantum discrete kinetic model [13-14] was performed and will be pre- 
sented here. In present approach the Uhling-Uhlenbeck collision term [13] which could describe 
the collision of a gas of dilute hard-sphere Fermi- or Bose-particles by tuning a parameter 7 : a 
Pauli-blocking factor (or 7/ with / being a normalized (continuous) distribution function giving 
the number of particles per cell) is adopted together with a disorder or free-orientation {6 which 
is related to the relative direction of scattering of particles w.r.t. to the normal of the propa- 
gating plane-wave front) into the quantum discrete kinetic model which can be used to obtain 
dispersion relations of (plane) sound waves propagating in quantum gases. We then study the 
quantum critical behavior based on the acoustical analog [15-17] which has been verified before. 
The possible phase diagram for MIT and/or resistance(or disorder)-scattering amplitude curves 
(as the temperature is decreased ) we obtained resemble qualitatively those proposed in [5] (cf 
Fig. 4 therein). The non-Fermi liquid behavior was also clearly illustrated here. 
We firstly introduce the concept of acoustical analog [15] in brief. In a mesoscopic system, where 
the sample size is smaller than the mean free path for an elastic scattering, it is satisfactory 
for a one-electron model to solve the time-independent Schrodinger equation : — (h 2 /2m)V 2 ip + 
V'(r)ip = Etp or (after dividing by —fi 2 /2m) V 2 4> + [q 2 — V{r)\^ = 0, where q is an (energy) 
eigenvalue parameter, which for the quantum-mechanic system is \j2mE/fi 2 . Meanwhile, the 
equation for classical (scalar) waves is V 2 ip — (d 2 "^ /(? dt 2 ) = or (after applying a Fourier 
transform in time and contriving a system where c (the wave speed) varies with position r) 
X7 2 ip + [q 2 — V(r)]ip = 0, here, the eigenvalue parameter q is uj/co, where to is a natural frequency 
and Co is a reference wave speed. Comparing the time dependencies one gets the quantum and 
classical relation E = fiuj [15-17]. 

We assume that the gas is composed of identical hard-sphere particles of the same mass [13, 14, 
18]. The velocities of these particles are restricted to, e.g., : ui, 112, ■ • • , u p , p is a finite positive 
integer. The discrete number densities of particles are denoted by Ni(x,t) associated with the 
velocities Uj at point x and time t. If only nonlinear binary collisions and the evolution of iVj 
are considered, we have 

^+u i -VN i = F i = ^(J%N k N l -A§N i N j ), i € A = {1, • ■ ■ ,p}, (1) 

j,k,l 

where (i,j) and (k,l) (i 7^ j or k 7^ I) are admissible sets of collisions [13-14,16-18] Here, the 
summation is taken over all j, k, I G A, where A 1 ^ are nonnegative constants satisfying [13-14,18] 
A J k \ = A l l x = AH, A ki( u i + Uj - u fc - uj) = 0, and A 1 ^ = Af y The conditions defined for the 
discrete velocities above require that there are elastic, binary collisions, such that momentum 
and energy are preserved, i.e., u« + uj = + u^, |uj| 2 + |uj| 2 = |u^| 2 + |u;| 2 , are possible for 
1 < k,l < p. We note that, the summation of N{ Ni) : the total discrete number density 
here is related to the macroscopic density : p (= m p J^i^i), where m p is the mass of the particle 
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[18]- 

Together with the introducing of the Uhling-Uhlenbeck collision term [13] : Fi = J2j,k,i ^ki WkNi 
(l+7iV i )(l +7JVj) -N i N j (l+'yN k )(l+'yN l )], into equation (1), for 7 < (normally, 7 = -1), 
we can then obtain a quantum discrete kinetic equation for a gas of Fermi-particles; while for 
7 > (normally, 7 = 1) we obtain one for a gas of Bose-particles, and for 7 = we recover the 
equation (1). 

Considering binary collisions only, from equation above, the model of quantum discrete kinetic 
equation for Fermi or Bose gases proposed in [13-14] is then a system of 2n(= p) semilinear 
partial differential equations of the hyperbolic type : 

a 2n 

at ax n 



d d cS 

7» Ni + v * " 7w N * = —T, N J N J+n(l + -yN j+1 )(l + 7^+n+l) 



2cS , JV i JV i+n (l+7JV i+ i)(l + 7JV i+n+ i), i = l,---,2n, (2) 

where iVj = Ni + 2 n are unknown functions, and v« =c(cos[# + (i — l)7r/n], sin[0 + — l)7r/n]); c is 
a reference velocity modulus and the same order of magnitude as that (c, the sound speed in the 
absence of scatters) used in Ref. 9, 9 is the orientation starting from the positive x— axis to the 
u\ direction and could be thought of as a parameter for introducing a disorder [15, 10, 16, 17, 20], 
S is an effective collision cross-section for the collision system. 

Since passage of the sound wave will cause a small departure from an equilibrium state and 
result in energy loss owing to internal friction and heat conduction, we linearize above equations 
around a uniform equilibrium state (particles' number density : Nq) by setting Ni(t,x) =Nq 
(1 + Pi(t, x)), where Pj is a small perturbation. After some similar manipulations as mentioned 
in [14-17], with B = 'jNq < [13, 14], which gives or defines the (proportional) contribution 
from the Fermi gases (if 7 < 0, e.g., 7 = — 1), we then have 

r d 2 2 2r„ (m-lk, a 2 „ „» T , n J ln 4cSN (l + B) " d „ 

fc=l 

where D m = (P m + P m+n )/2, m = 1, • • • , n, since Di = D m for 1 = m (mod 2n). 
We are ready to look for the solutions in the form of plane wave D m = a m exp i(kx — ut), 
(m = 1, ■ ■ ■ ,n), with u;=a;(fc). This is related to the dispersion relations of ID (forced) plane 
wave propagation in Fermi gases. So we have 

(1 + ih (i + B )- 2X 2 cos 2 [9 + (m ~ 1)?r ])a m - V flfe = 0, m = 1, • • • , n, (4) 

n n ^ 

where 

A = kc/(V2uj), h = AcSN /uj oc 1/K„, 

where /i is the rarefaction parameter of the gas; K n is the Knudsen number which is defined as 
the ratio of the mean free path of gases to the wave length of the plane (sound) wave. 
We can obtain the complex spectra (A = A r + i A^; A r = k r c/ (\f2~uj): sound dispersion, a relative 
measure of the sound or phase speed; \ = fcjc/ (V2lo) : sound attenuation or absorption) from 
the complex polynomial equation above. Here, B could be related to the occupation number of 
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different-statistic particles of gases To examine the critical region possibly tuned by the Pauli- 
blocking measure B = jNq and the disorder 6, as evidenced from our preliminary results : 
Xi = for cases of B = — 1 or 9 = it/A [17], we firstly check those spectra near 9 = 0, say, 
9 = 0.005 and 9 = tt/4 ~ 0.7854, say, 9 = 0.78535 for a B-sweep (B decreases from 1 to -1). We 
plot them into figures 1, 2, respectively. Note that, as the disorder or free-orientation 9 is not 
zero, there will be two kinds of propagation of the disturbance wave : sound and diffusion modes 
[16-17,19-20]. The latter (anomalous) mode has been reported in Boltzmann gases [19,21] and 
is related to the propagation of entropy wave which is not used in the acoustical analog here. 
The absence of (further) diffusion (or maximum absorption) for the sound mode at certain state 
(h, corresponding to the inverse of energy E; cf. Refs. 11 or 15) is classified as a localized state 
[11,15,17]. The state of decreasing h corresponds to that of T (absolute temperature) decreasing 
as the mean free path is increasing (density or pressure decreasing). 

We can observe the max. Aj (absorption of sound mode, relevant to the localization length 
according to the acoustical analog [15,17]) drop to around four orders of magnitude from 9 = 
0.005 to 0.78535! This is a clear demonstration of the effect of disorder. Meanwhile, once 
the Pauli-blocking measure (B) increases or decreases from zero (Boltzmann gases), the latter 
(Fermi gases : B < 0) shows opposite trend compared to that of the former (Bose gases : B > 0) 
considering the shift of the max. Aj state (5h). 5h > is for Fermi gases (\B\ increasing), and 
the reverse (Sh < 0) is for Bose gases (B increasing)! This illustrates partly the electron-electron 
interaction effect (through the Pauli exclusion principle) . These results will be crucial for further 
obtaining the phase diagram (flow to metallic or insulating state as the density or temperature is 
decreased) tuned by both disorder and the (electron-electron) interaction below. Here, B = — 1 
or 9 = 0,7r/4 might be fixed points commented in [22]. 

To check what happens when the temperature is decreased (or the density or h is decreased) to 
near T = 0, we collect all the data based on the acoustical analog from the dispersion relations 
(especially the absorption of sound mode) we calculated for ranges in different degrees of disorder 
(here, is up to 7r/4 considering single-particle scattering and binary collisions; in fact, effects 
of 9 are symmetric w.r.t. 9 = it/A for < 9 < 7r/2 [17]) and Pauli-blocking measure. After 
that, we plot the possible phase diagram for the (dimensionless) conductivity vs. the absolute 
(dimensionless) temperature into Fig. 3 (for different Bs : B = —0.9,-0.7,-0.5,-0.3,0.1). 
Here, MFP is the mean free path and the temperature vs. MFP relations could be traced from 
[23] (cf. Fig. 3 therein). Note that, the resistivity or resistance is proportional to the strength of 
disorder in 2D (in the sense that for weak disorder it is given by 1/(&fZ), in units of h/e 2 ; kp is 
the Fermi wave number, I is the mean free path associated with the usual Drude conductivity) 
[5]. This figure shows that as the temperature decreases to a rather low value, the resistivity (or 
strength of disorder) will decrease sharply (at least for Bose or Fermi gases). There is no doubt 
that this result resembles qualitatively that proposed before [5,24-25]. 

To know the detailed effects of electron-electron interactions (tuned by Bs here), we plot the 
corresponding figure as shown in Fig. 4. Each contour line (flow path relevant to a specific B 
tuned by the disorder) represents the behavior when T is decreased, and different (flow) lines 
represent different values of scattering amplitude (Sb oc K n or r s ). Interesting results are (i) 
there seems to be a quantum phase transition boundary (interface or regime) for bosonic-like 
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particles (B > 0) which resembles that of high-temperature superconducting phase transition 
due to doping if we treat the strength of disorder to be equivalent to the doping amount! (ii) 
as evidenced in the top part of this plot for the Bose gases (near QCPs for smaller resistivity 
together with rather large arrows), it confirms Larkin's comment (cf page 793 in [2] by Larkin) : 
both Bose and Fermi approached are important for the (QP) transition; Bose approach is more 
useful in a small region close to the transition! 

To make sure we already recover previous proposed results (possible renormalization group (RG) 
flows for disorder plus interactions or scattering amplitude-resistance curves or suggested phase 
diagram tuned by r s and disorder, cf. [5] or [24-25]), we summarize our results by illustrating 
them into a 3D plot as demonstrated in Fig. 5. Crossover lines separate those flows which begin 
at T > at small resistance (disorder or the inverse of conductance : 1/G) and flow initially, 
as the temperature is lowered, toward larger resistance but then are repelled by the fixed point 
and flow to large Sf, (or r s ) and large G (metallic behavior) from those which begin at larger 
disorder (or smaller G) and flow toward very-large disorder (insulating behavior). Again, this 
result resembles that proposed before (cf. [5] or [24-25] therein). Note that, this flow is also 
similar to those in [24-25] (e.g., cf. Fig. 41 by Aoki [24]). Meanwhile, as expected before, at 
T = 0, a 2D system would become a Wigner crystal phase at r s > 37 [26]. This lies possibly 
near or above the rather large Sb position with zero strength of disorder in our illustrations. 
Possible QCPs happen around S b ~ 80 for B = 0.85 but S h ~ 10 for B = 0.1. 
To conclude in brief, our illustrations here, although are based on the acoustical analog of 
our quantum discrete kinetic calculations, can indeed show the non- Fermi liquid and quantum 
critical behavior for the metal-insulator transition in 2D [27-31] as the temperature is decreased 
to rather low values. 

References 

[1] J. Custers, et al, Nature 424, 524 (2003). J. Hertz, Phys. Rev. B 14, 1165 (1976). A.J. 
Millis, Physica B 312-313, 1 (2002). 

[2] A.M. Goldman, Physica E 18, 1 (2003). A. Larkin, Ann. Phys. (Leipzig) 8, 785 (1999). 
M.P.A. Fisher, Phys. Rev. Lett. 65, 923 (1990). 

[3] M. Greiner, et al, Nature 415, 39 (2002). 

[4] G. Modugno, et al, Science 297, 2240 (2002). K.E. Strecker, G.B. Partridge, and R.G. 
Hulet, Phys. Rev. Lett. 91, 080406 (2003). C.A. Regal and D.S. Jin, Phys. Rev. Lett. 90, 
230404 (2003). 

[5] E. Abrahams, Physica E 3, 69 (1998). 

[6] E. Abrahams, S.V. Kravchenko, and M.P. Sarachik, Rev. Mod. Phys. 73, 251 (2001). 

[7] L.D. Landau, Sov. Phys. JETP 5, 101 (1957). CM. Varma, Z. Nussinov, and W. van 
Saarloos, Phys. Rep. 361, 267 (2002). 

[8] N.F. Mott, Rep. Prog. Phys. 47, 909 (1984). P.W. Anderson, Comments Solid Stat. Phys. 
2, 193 (1970). 

[9] A.M. Finkel'stein, Z. Physik B 56, 189 (1984). B.L. Altshuler, D.L. Maslov, and V.M. 
Puladov, Physica E 9, 209 (2001). A.P. Taylor and A. MacKinnon, J. Phys. Condens. 
Matter 14, 8663 (2002). V. Dobrosavljevic, et al, Phys. Rev. Lett. 79, 455 (1997). C. 
Castellani, et al, Phys. Rev. B 30, 1596 (1984). 



6 



[10] P.W. Anderson, Phys. Rev. 109 1492 (1958). A. Figotin and A. Klein, Commmun. Math. 
Phys. 180 439 (1996). 

[11] T.R. Kirkpatrick, Phys. Rev. B 31 5746 (1985). 

[12] S. John, Phys. Today 44 32 (1991). A. Klein and A. Koines, Mathematical Physics, Analysis 
and Geometry 4 97 (2001). M. Aizenman, J.H. Schenker, R.M. Friedrich, and D. Hundert- 
mark, Commun. Math. Phys. 224 219 (2001). 

[13] V.V. Vedenyapin, I.V. Mingalev, and O.V. Mingalev, Russian Academy of Sciences Sbornik 
Mathematics 80 271 (1995). 

[14] A. K.-H. Chu, J. Phys. B At. Mol. Opt. Phys. 34 L711 (2001). A. K.-H. Chu, Phy. Scripta 
(to be appeared, 2004). 

[15] J.D. Maynard, Rev. Mod. Phys. 73 401 (2001). 

[16] K.-H. W. Chu, J. Phys. A Math, and General 34, L673 (2001). W. K.-H. Chu, Appl. Math. 
Lett. 14 275 (2001). 

[17] A. K.-H. Chu, Phys. Rev. E 66 047106 (2002). 

[18] T. Platkowski and R. Illner, SIAM Rev. 30 213 (1988). N. Bellomo and T. Gustafsson, Rev. 
Math. Phys. 3 137 (1991). 

[19] A. K.-H. Chu, Eur. Phys. J. B 10 1 (1999). 

[20] A. K.-H. Chu, Preprint (2002). A. K.-H. Chu, (in Press, 2003). 

[21] R.J. Jr. Mason, in Rarefied Gas Dynamics, edited by J.H. de Leeuw (Academic Press, New 
York, 1965), vol. 1, p. 48. 

[22] P.W. Anderson, Physica B 318, 28 (2002). 

[23] D. Einzel and J.M. Parpia, J. Low Temp. Phys. 109 1 (1997). 

[24] H. Aoki, Rep. Prog. Phys. 50, 655 (1987). 

[25] D.V. Shopova and D.I. Uzunov, Phys. Rep. 379, 1 (2003). 

[26] E. Wigner, Phys. Rev. 46, 1002 (1934). B. Tanatar and D.M. Ceperley, Phys. Rev. B 39, 
5005 (1989). 

[27] A. Punnoose and A.M. Finkel'stein, Phys. Rev. Lett. 88, 016802 (2002). M. Hilke, Phys. 
Rev. Lett. 91, 226403 (2003). 

[28] V.M. Pudalov, et al, Phys. Rev. Lett. 91, 126403 (2003). A. A. Shashkin, et al, Phys. Rev. 
Lett. 86, 086801 (2001). 

[29] Y.Y. Proskuryakov, et al, Phys. Rev. Lett. 89, 076406 (2002). R.J. Bell, Rep. Prog. Phys. 
35, 1315 (1972). 

[30] R. Shankar, Rev. Mod. Phys. 66, 199 (1994). S. Sachdev, Quantum Phase Transitions 
(Cambridge Univ. Press, 1999). T. Chakraborty and V.M. Apalkov, Adv. Phys. 52, 455 
(2003). A.F. Hebard and M.A. Paalanen, Phys. Rev. Lett. 65, 927 (1990). 

[31] M.P.A. Fisher and D.H. Lee, Phys. Rev. B 39, 2756 (1989). 



7 




10"' 10° 10 1 



Rarefaction measure : h « 1/K 
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Fig. 1 Variations of the smaller (Sound mode) and larger (Diffusion mode [19,21]) 
A, w.r.t. h for the same disorder : 9 = 0.005 in different Pauli-blocking measures 
(B > corresponds to Bose gases). The maximum absorption states for sound modes 
correspond to possible localized states as the diffusion stops. B decreases from 1 

(the highest) to —1 with the decrement being —0.2. Xi = as B = —1. h = 4cSN /lo. 

N is the number density and S is the effective collision cross-section. 




B= - 0.8 



Rarefaction measure : h « 1/K 

n 

Fig. 2 Pauli-blocking (B) effects on the absorption Xi. (Note that the 
localization length oc 1/A i; cf. [11,16-17]. The energy E (relevant to the 
illustration of localized states in [11]) corresponds to Tilo [15-16]. E oc 1/h 
once cS No is fixed.) B decreases from 1 (the highest) to —1 with the 
decrement being —0.2. Bose and Fermi gases are in opposite trend for 
the same increment \B\ = 0.2. 6 = 0.78535 and A, = as 6 = tt/4 k, 0.7854. 
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Fig. 3 Possible phase diagram for non-Fermi gases w.r.t. the disorder and 
the scattering amplitude £&. <S& oc BK n and if n oc the mean free path. Note that 
the resistivity is proportional to the strength of disorder in 2D [5-6]. Here, once 
the temperature is lowered down the resistivity firstly increases and then decreases. 
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Flow (arrow) direction : Temperature decreases 
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Fig. 4 Possible flow pattern for Sj, w.r.t. the disorder and Pauli-blocking 
measure (B). Sb oc BK n and K n oc the mean free path. Note that the 
resistivity is proportional to the disorder [5-6] in 2D. Each continuous flow 
represents the behavior when T is decreased, and different paths represent 
different values of Pauli-blocking measure, (cf Fig. 4 in [5]) 




Fig. 5 Possible scattering amplitude (St,) vs. disorder and Pauli-blocking 
measure curves. Sj, oc BK n , and K n oc the mean free path. Note that 
the resistivity is proportional to the strength of disorder in 2D [6]. Here, 
the flow (arrow) direction means that the temperature is decreasing. 
These results resemble qualitatively those proposed by Abrahams in [5-6]. 



